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Abstract 

The radiation emission spectra of polarized photons emitted from charge particle in 
single crystal are obtained in semiclassical theory in Baer-Katkov-Strakhovenko approx- 
imation for planar crystal orientation. The range of applicability of this approximation 
is estimated by comparing the results with calculations in exact semiclassical theory. 
Optimal crystal orientations for producing unpolarized and polarized photon beams are 
also founded. 

1. Introduction. It is well known that the processes of photon emission from charged 
particles (as well as e~e + pair creation by photon) at high energies in oriented single crystals 
reveal very strong angular and energy dependence. The photon emission spectra in crystals 
in their maxima could exceed the corresponding values for amorphous target 100 and more 
times and the radiation appears essentially polarized. Therefore, these processes are widely 
applied in experimental physics for the production of high energy polarized photon beams, as 
well as for the analyses of photon polarization. For incident angles •& to crystal axis/planes 
i? $ v , where $ v is the characteristic angle given by $ v = U /m (U is the scale of ax- 
ial/planar potential and m the electron mass) these processes are well described in the theory 
of coherent bremsstrahlung (CB) and pair production (CPP) [1]. This theory is constructed 
in the framework of the first Born approximation in crystal potential and fails at small angles, 
where these processes become of magnetic bremsstrahlung nature. 

In the papers of V.N. Bayer V.M. Katkov and V.M. Strakhovenko has been developed 
a general theory of basic QED processes in strong crystalline fields (see i.e. new English 
edition of their book [2]). This theory is not restricted by the first Born approximation and 
is based on the semiclassical character of motion of ultrarelativistic particles in strong fields. 
However, the numerical calculations in the exact theory face serious mathematical difficulties 
and some approximation and modeling methods have been developed. First numerical results 
of the radiation emission spectra in the semiclassical theory were published in [3] (in the 
framework of additional modeling assumption); exact calculations were presented in [4]. In 
Ref. [5] an analytical method of calculations in the framework of semiclassical approach 
was developed and new numerical results were presented. The radiation emission spectra of 
polarized photons in semiclassical approach were calculated in [2,6] and first numerical results 
were presented in our previous paper [6]. These papers demonstrate an essential difference 
of exact spectra from the corresponding results of coherent theory at small angles $o < 
There is a good agreement between first experimental results at small angles [7,8] and these 
calculations. 

Polarization characteristics of radiation emission in single crystals in the region of appli- 
cability of CB are well investigated in experiments (at the energies of incident electrons up 
to 10 GeV) and the CB theory describes the experimental results satisfactory [1]. 

1 The work is supported by ISTC grant A-099 
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2. Main results. In this paper we have derived the polarized photon emission spectra 
in the framework of Bayer-Katkov-Strakhovenko approximation of paper [3] (further BKS 
approximation) which is valid at small incident electrons angles to one of crystallographic 
planes. We also modified the nonpolarized equations of [3] and present all results in the form 
very convenient for numerical calculations. Our calculation algorithm allows to reduce the 
computing time of all polarized spectra to 5-10 sec (Pentium 233) depending on accuracy and 
crystal orientation. 

Comparing the calculation in BKS approximation [4] with exact calculations in semiclassi- 
cal theory we estimated the range of applicability of this approximation. It appears, that the 
equation in BKS approximation described the exact spectra satisfactory only at very small 
angles of incidence, about 3-5 times less than critical angle of planar channeling. 

We present a large number of polarized emission spectra in wide range of angles for 
different energies and crystal orientations. The analysis of these results make it possible to 
find optimal crystal orientations for producing high energy photon beams. For producing 
nonpolarized high energy photons it seems optimal to use the range of hard photon peaks at 
large x, where x is a fraction of electron energy e carried by emitted photon. The maximum 
of intensity we find for < 110 > oriented crystal for incident electron moving close to (111) or 
(110) planes. However, the degree of polarization at hard photon peaks appeared to be rather 
low. We find that the degree of polarization at given x achieves its possible maximal value 
(~ 55% at x—2/3, at e — 200 GeV) at relatively large values of incident angles to crystal 
axes, when the corresponding hard peak situated at given x is of pure coherent nature. The 
optimal configurations for producing polarized high energy photons we find for < 001 > 
oriented diamond crystal when electron moves close to (110) plane and for < 110 > oriented 
crystal when electron also moves close to (110) plane. 

3. Polarized photon spectrum in BKS approximation. Let us start from the 
general semiclassical formula for polarized photon spectrum [2.6], which we present in the 
following form in terms of Stocks parameters £i, £ 2 , 
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Here e and u are electron and photon energies, 7 = e/m; r is the entry point, V is crystal 
volume, integration over r gives the averaging over electron trajectories in crystal, F(r, 1?) is 
the coordinate distribution function at given value of Eq.l must be also averaged over time 
t. The quantities of Ai = Ai(t) and 7] = Tj(r),i = 0, 1,2,3, in (1) in rectilinear trajectory 
approximation have the following form: 
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q y q' x )[g(q z r) sm(q' z r) - g(q' z r) sm(q z r]e-^^\ (5) 

q y q' y )[g{q z T)g{q' z T) + sm(q z r) Bmtf z T)]e-**«™ , (6) 
sin a; 

q(x) = cosx. 

The notation used here coincides, in general, with ones of Ref . [2] . The crystal potential is 
used in the form of U(r) = J2G(q)e~ iqr , where q is reciprocal lattice vector; e' = e — u, 
tp{e) — j + jf, r(t) = r + v t The Stocks parameters and projections of vectors correspond 
to the frame defined by following unit vectors: 

e x = [n 2 v ], e y = n 2 , e z = v , (7) 

where vo is the incoming electron velocity, ri2 is a unit vector transverse to vo and to the 
main crystallographic axis. 

Let d be the electron incident angle to the given crystallographic axes and ip be electron 
incident angle to the given crystallographic plane, then sin^ = sin d sin tp, where tp is the 
angle of v projection onto the plane perpendicular to the given crystallographic axis, to the 
given crystallographic plane. 

The BKS approximation is valid for planar crystal orientations under the following con- 
ditions. 

1. tp <C 1, this makes it possible to isolate the continuous planar potential contribution 
splitting the sums in (2-6) into two parts, J2q — Sq + Dq\ where the sum J2q contains 
q for which q z — > when tp — > and the sum J2q includes all other terms. 

2. ^<l,^i should be small enough, so that \q z \r <^ 1, this makes it possible to expand 
the functions in the sum J2q m powers of \q z \r. 

3. $ should not be too small, ■& ^> i?„ so that the sum J2q be small enough and Eqs.2-6 
could be expanded in powers of ^q- 

Two last conditions crucially depend on the effective upper integration limit over r in (1), 
To. The value of tq could not be exactly established from general considerations due to rather 
complicated dependence of Eq.2 on r, so, further we try to estimate the range of validity of 
BKS approximation by comparing the BKS results with calculations in exact semiclassical 
theory. 

Following the steps of Ref. 3, we did the above mentioned expansions, performed time 
averaging and integrated the results over r. All polarized intensity spectra we succeeded to 
present via one fold integrals of known special functions: 
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where the integration is taking over coordinate y perpendicular to the plane, d p i is the distance 
between the planes. Circular polarization vanishes in BKS approximation, I 2 — 0. The 
quantities I , I±, J 3 are defined as follows: 
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u = u/e', s = 2eq z /m 2 ) A 2 = x 2 «, (3 = f{e), z = u(l + p/2), z± = u(l + p/2 ± s/u). 
The functions F (z, A), F%(z, A) and F 2 (z, A) could be expressed via known special functions. 
We present and analyze them in Appendix. The nonpolarized intensity spectra, Jo, coincide 
with result of Ref.[2], if we take for F (z,X) the integral representation via Bessel functions 
(Eq.7of[2]). 



c<1 1 0>, in (1-11) plane 
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4. Numerical calculation results and conclu- 
sions. In this section we present the numerical calcula- 
tions for uniform distributions, keeping, of course, in 
mind, that for real spectrum one should take into ac- 
count the particle flux redistribution (F(y, 9q) 7^ 1) for 
ip < ip c , where ipc is the critical angle of planar chan- 
neling, if) c = yf^f, Uq is the depth of planar potential 
well. For < 001 > orientation we define the angle <p as 
the angle of the projection of v onto the plane (001) 
to the plane (110), similarly, for < 110 > orientation 
the angle <p is the angle of the projection of v onto 
the plane (110) to the plane (001). The typical inten- 
sity spectrum and Stocks parameters are presented on 
Fig.l for < 110 > oriented diamond crystal. Note, that 
throughout the paper we present all angles in Radian. 
The soft part of the spectrum is completely determined 
by magnetic bremstrahlung contribution. There is a 
wide polarization maximum of essential linear polariza- 
tion (£i = 0, £3 ~ 0.5 — 0.7) in wide range of x < 0.3. 
The shape of intensity and polarization spectra in soft 
region are not changed by the variation of the angle (p, 
because in BKS approximation the continuous planar 
potential contribution was actually derived in the limit 
if —>■ 0. The high frequency peaks in hard part of the 
spectra could be interpreted as coherent peaks slightly 
modified by influence of continuous planar potential. 

This region can be approximately described by modified coherent theory [2], if we exclude 
from the sums the reciprocal vectors corresponding to magnetic bremstrahlung contribution 
(q z — > when ip — > 0). It is demonstrated on Fig.l. In such approximation each hard peak 
corresponds to the sum of coherent contributions of reciprocal vectors, so that for each peak 
corresponding values of q z tend to the fixed limit at ip — > [6]. The polarization in hard 
region strongly depend on crystal orientation and will be discussed later. 

In order to estimate the range of applicability of BKS approximation we also perform ex- 
act semiclassical calculations (making use of Eq.1-6) at different incident angles to the plane 




X =co/e 



Fig.l: Comparison the hard part of 
spectral distributions of intensity and 
Stocks parameters in modified coherent 
and semiclassical (BKS) approaches 
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C<001>, E = 200GeV, e=0.0002 



0.02 (BKS approx.) 
0.075 (chan.crit. angle) 
0.04 
0.02 





C<001>, E=600GeV, e=0.0002 



<p=0.02 (BKS approx.) 
<P=0.075 

<P=0.043 (chan. crit. angle) 
<P=0.02 





X = (o/E 



X=co/E 



Fig. 2: Comparison the spectral distributions of intensity and degree of polarization of BKS approx- 
imation with exact semiclassical calculations. 

for different energies. To carry out the arising integrals over time of complicated oscillating 
function we have elaborated special integration program. Results are presented in Fig.2. It 
is very interesting that when electron incident angle to the plane is equal or even less that 
planar channeling critical angle, ij) c , the soft part of the spectrum strongly differ from BKS 



C<001>, cp=0.322, in (-130) plane 



0-0.00046 , E.150 GeV 
0.0.00025, E.300 GeV 
0.0.00090, E= 80 Gev 




Fig.3: Energy dependence of spectral distributions of intensity and degree of polarization. 

calculations and is far from pure magnetic bremstrahlung nature. It more resembles the 
modified coherent spectrum, with typical peaks and dips. That means that at these values 
of incident angles the functions in the sums J2q cannot be expanded in powers of \q z \r in 
all integration region, and together with magnetic bremstrahlung contribution, which comes 
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from low t region also contribute coherent effects which come from high r region. Only at 
very small angles of incidence, about 3-5 times ( for e ~ 100 — 600GeV) less than if) c , the soft 
part of spectrum become of pure magnetic bremstrahlung nature and could be described in 
BKS approximation. Note, that hard part of spectra is well described in BKS approximation 
starting from large values of incident angles (practically from $ < 
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Fig. 4: Spectral distributions of intensity and degree of polarization for different crystal orientations 
for fixed position of hard photon peak. 

We also present large number of intensity and polarization spectra to demonstrate the de- 
pendence of angles, energies and crystal orientations and make it possible to find optimal 
configurations for producing high energy photon beams. Some of the analogous nonpolar- 
ized intensity spectra were already presented and analyzed in different papers [2-6]. The 
energy dependence at fixed incident angles is demonstrated in Fig. 3a. Hard peaks in inten- 
sity spectra move to the left with increasing energy, slightly change their form, the heights 
of the peaks being practically independent of energy, corresponding peaks in polarization 
spectra also moved to the left and their heights decrease. Soft peak also moves to the left 
and becomes wider. On Fig. 3b we demonstrate the energy dependence of spectra at different 
angles, chosen in such a way that first hard peak is situated at fixed value of x ~ 2/3. The 
height of hard peak in intensity spectra increase with increasing energy, the hard part of 
polarization spectra practically remain unchanged. So at very high energies it will be more 
and more preferable to use high frequency peaks in photon spectra for producing high energy 
photon beams. On Fig.4 we present the spectra for < 001 > and < 110 > diamond crystal 
orientation when electron moves close to the main crystallographic planes at fixed energy. 
Angles are arranged in such a way that first hard peak in all spectra is situated at x ~ 2/3. 
The maximal height of intensity peak takes place for < 110 > orientation when electron moves 
close to (111) or (110) planes and these configurations are the most preferable for production 
of nonpolarized high energy photons. The degree of polarization in hard photon peaks for 
these configurations appeared rather low. It is interesting that if we are going from strong 
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planes to weak ones, the intensity of hard photon peak is decreasing and degree of polar- 
ization increasing. All these features could be qualitatively understood in modified coherent 
theory. At hard photon peaks the contributions of different planes are summed in nonpo- 
larized photon intensity (J ), but in polarized photon intensities (Ji, J 3 ) they could have 
opposite signs or one of the planes could not contribute. Therefore, the degree of polarization 
very strongly depend on crystal orientation. As one can see from Fig. 4, for some orienta- 
tions different contributions in polarized intensities practically compensate each other. It 
seems impossible to find crystal orientation, for which the contributions of all main plains 
will be summed in I-y and I3 . Therefore, the maximal polarization for hard photon peak 
at given value of x can be achieved for sufficiently large angles of incident electron to the 
crystallographic axis when hard photon peak appeared to be a typical coherent peak, corre- 
sponding to contribution of one plane and isn't affected by influence of even very weak planes. 
This is also demonstrated in Fig. 5 We present the 
intensity and polarization spectra for < 001 > dia- 
mond crystal for electron moving close to (110) plane 
for large 9 and different <p, keeping the position of 
main peak fixed at x ~ 2/3. In wide range of (p , 
starting from p ~ 7r/4, height of the peak and degree 
of polarization remain practically unchanged and co- 
incide with presented curves for p =7r/4-0.05, if we 
keep electron far from neighboring planes. When we 
come closer to even rather weak planes, the height 
of the peak increases and the degree of polarization 
decreases. We demonstrate this effect for weak (350) 
plane (Uq ~ 1.2 eV). It is worth also noting that when 
we come close enough to such week plane, we could see 
even considerable magnetic bremstrahlung contribu- 
tion. We have also checked this properties for weaker 
planes, (370) , (570), (390) and have found that their 
influence cannot be entirely ignored. 

The most preferable crystal orientations for pro- 
duction of polarized high energy photons we find for Fig.5: The influence of high order plane 
< 001 > oriented diamond crystal when electron on intensity and polarization in coherent 
moves close to (110) plane (for instance at 9 = 0.028 peak, 
rad, ip = 0.735 rad for x ~ 2/3) and for < 110 > 

oriented crystal when electron also moves close to (110) plane (for instance, at 9 = 0.028 rad , 
<p = 1.536 rad for x ~ 2/3). For these configurations the maximal degree of polarization could 
be achieved for maximal intensities. These configurations bring exactly to the same results. 
The exact values of angles do not matter, one only should arrange angles to fix the main 
coherent peak at given value of x and keep angles far from influence of neighboring planes. 

This work was supported by the ISTC grant, project A-099. One of us (N.T) is grateful 
to Prof. R.Avakian for suggestion to participate in this project. 
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Appendix 



The function F (z, A) from Eqs.(9-ll) can be expressed through generalized hypergeomet- 
ric functions, p F q ((a p ); (b p ); x): 

Fo(Z ' A) "6 + (3A)W(|) 4A^ (A ' 1} 

Functions A) and F 2 (z, A) are expressed through combinations of Bessel functions of 
the first kind, J v {x) , and of the third kind, K v (x): 
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Here is Euler gamma function and is Heaviside function. Note, that: 

Fl(t , A) = ^), F2(2 , A) = _^>. 

Instead of making use of existing programs for calculations of special functions, we have 
elaborated the calculation algorithm directly for functions F^z, A), connecting the asymptotic 
expansion in powers of \ jz with low z expansion in power of z. We present the corresponding 
equations only for F (z,X), analogous equations for F 2 (z,X) and Fi(z,X) could be derived 
by differentiation or from corresponding expansions of Bessel functions. The asymptotic 
expansion of F (z, A) has the following form: 
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Power series expansion in powers of z follows from definition of generalized hypergeometric 
functions: 
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where t = z 3 /3X 2 . Asymptotic expansions (A.4-A.5) work with accuracy better than 0.1% 
starting from \z\ /A 2 / 3 ~ 4 and Eqs. (A.4-A.5) and (A.6) could be connected at this point. 
One should take at this point about 10 terms of expansion (A.6) . Such calculation algorithm 
proved to be very effective and essentially reduces the computer calculation time. 
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